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We discuss some problems concerning the application of pcrturbative QCD to high energy 
processes. In particular for hard processes, we analyze higher order and higher twist 
corrections. It is argued that these effects are of great importance for understanding the 
behaviour of pion electromagnetic form factor at moderately large momentum transfers. 
For soft processes, we show that summing the contributions of the lowest twist operators 
leads to a Regge-like amplitude. 
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1. Introduction. 

It is now safe to assert that quantum chromodynamics (QCD) agrees qualitatively 
with all the experimental data related to strong interactions phenomena. At the 
same time, the QCD predictions are usually too flexible for precisive quantitative 
tests. This is caused mainly by the fact that all the calculations in QCD are based on 
perturbation theory (PT) , i.e. on the expansion over the coupling "constant" a s (k) 
that depends really on the momentum scale k related to the process investigated. 
Asymptotic freedom^ enables one to use PT at short distances (or large momenta) . 
However, any physical process involves also long distances, i.e. each process involves 
small momentum scales p 2 (e.g. , quark and hadron masses) , and as a rule, this results 
in the appearance of the logarithmic contributions hi(Q 2 /p 2 ), that are singular for 
p 2 = (mass singularities^™^) . In such a situation » 2 cannot be neglected. However, 
within PT it ispossible to show that for inclusive^-™^ and some hadronic exclusive 
hard processea^™ the Q 2 -dependence of the corresponding amplitude T(Q 2 ,p 2 ) 
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mitted to the XX International Conference on High Energy Physics, Madison, July 1980. 



1 



November 6, 2009 8:3 Efremov 




2 A. Efremov and A.Radyushkin 

can be factorizecjf] from the p 2 -dependence (see Fig|TJ): 

T(Q 2 ,p 2 ) = Q N {E(Q 2 /fi 2 , a s (fi)) f(fi 2 , P 2 ) + R(Q,p)} , (1.1) 

where N is the dimension of T in mass units and R is the sum of contributions 
which are power suppressed with respect to E <8> / ■ The parameter [i is a boundary 
between large and small momenta, and E®f does not depend on a particular choice 
of fJL. 

The functions / describe 
long-distance interactions. This 
means they cannot be reli- 
ably calculated in perturbative 
QCD and must be treated phe- 
nomcnologically. As a result, 
the QCD predictions are more 
ambiguous. The functions E 
describe short-distance interac- 
tions. In principle, they are Fig . L Illustration of Equation CEQ. 
given by a perturbative series 

expansion over a s (/x). In practice, only a few terms are known, and one is forced to 
make some plausible hypotheses about the magnitude of the non-calculated higher 
order corrections. According to most of recent estimates, a s {/j,)/% is of the order 
0.1 for [i 2 < lOGcV 2 . This means that taking into account only a few first terms 
of the series is a good approximation only if the coefficients of the expansion of E 
over a s /ir are not too large. 

An analogous uncertainty exists also for power corrections absorbed by R(Q,p). 
It is known that the hadron and quark mass corrections for the most simple inc lusive 
processes can be calculated exactly with the help of the ^-scaling formalis m^ 
For light quarks u,d,s, their masses (m u ~ 4McV, rtid — 7McV, m s ~ 120 MeV) 
usually may be neglected at all. The main uncertainty is due to power corrections 
caused by the finite size of the hadrons, by the Fermi-motion of quarks inside the 
hadrons, etc. All these effects have a nonperturbative origin. The magnitude of 
the corresponding corrections (M 2 /Q 2 ) n is determined by a characteristic scale 
M ~ 1 /i?conf — 300 -T- 500 MeV, but in some cases they play a very important role 
up to very high momentum transfers. For example, in high-py hadron production, 
the effects of th e primordial transverse momentum of partons dominate the cross 
section up ■p 2 T ~ 30 -r- 40GeV 2 . It should be emphasized, however, that all 
the phenomenological methods of taking into account the power suppressed terms 
have no reliable theoretical basis. So, it is highly important to develop methods of 
field-theoretical analysis of power corrections. 



^{Comm. 2009.) For the Yukava theories (scalar gluons) the property of factorization for hard 
processes was first considered even earlier (A. V. Efremov, Yad. Fiz. 19, 196 (1974)). 
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2. Higher-order corrections 

In general, the functions E(Q 2 /fi 2 , a s (ft)) in Eq. (jl.lj) depend on the calculation 
scheme, namely, on the chosen recipe of the R-operation for the ultraviolet diver- 
gences and the recipe of separating the contributions related to short and long 
distances (i. e. on the R-operation for composite operators). In particular, E is 
/^-dependent. The functions f(fJ?, p 2 ) also depend on the chosen scheme and only 
the product E ® f is scheme-independent. If we take fi = Q, then the resultant 
expression would not have an explicit dependence on ft. Note, that this procedure 
removes from E{Q 2 / ^t 2 ) the logarithms ln(Q 2 /^ 2 ) which, for Q ^> ft,, lead to growth 
of the coefficients in the expansion of E over a s . The meaning of the choice (i = Q 
is clear: one must take ft equal to a scale characterizing the off-shellness of the par- 
ticles taking part in the short-distance subprocess and the latter is proportional to 
Q 2 : (k 2 ) ~ ~a 2 Q 2 . 

If, however, the parameter a is very large (or very small) compared to 1, then 
the choice ft ~ aQ should be preferred. It is implicit here that we use a "physical" 
renormalization scheme, i.e. g(k) corresponds to a vertex with external momenta k 2 . 
However, for direct calculations within QCD it is more convenient to use various 
"unphysical" schemes based on dimensional regularization. In this case the meaning 
of g(k) is less transparent. So, for the time being, we will adhere to the choice /i = Q. 

Recently, in a series of papers^HISl it has been established that the scheme- 
dependence of the results obtained may be reduced by expanding the coupling 
constant a s (Q) over (lnQ 2 /A 2 ) -1 , 

4 + (") 



4tt b Q \nQ 2 /A 2 \ b 2 lnQ 2 /A 2 

where bo, £>i, . . . are the coefficients of the expansion of the /3-function over g. After 
this change we have the following representation for E <£> /: 



(\nQ 2 /A 2 y>°/ bo 



(&olnQ 2 /A 2 )<, 

a 10 , a 20 



+ 



6 lnQ 2 /A 2 6 2 ln 2 Q 2 /A 2 
q 2 i lnlnQ 2 /A 2 
b 2 In 2 Q 2 /A 2 + ' 



(2.2) 



All information about the long-distance dynamics is accumulated in /, whereas 
the coefficients aik can be calculated in PT. Moreover, the coefficients , a^} 
related to two different schemes can be obtained from one another by the change 
Ai = K12A2 for the appropriately chosen parameter k,\2- Thus, if one uses the 
expansion ()2.2|) . then various schemes differ only in magnitude of the parameter A. 

Let us assume that, in a scheme S, some first coefficients aik are numbers of 
the order 1. Then in another scheme S", which has A' = 100 A ( or A' = 0.01 A), 
the coefficients aik are numbers of the order bo In 100 ~ 40. It is easy to notice the 
analogy with our previous discussion about the optimal choice of the parameter /i 



November 6, 2009 8:3 



Efremov 



4 A. Efremov and A.Radyushkin 

and to conclude that the scheme S is very close to a "physical" scheme, since the 
choice fi = Q (assumed in Eq. I|2.2p ) minimizes for this scheme the higher-order 
corrections. Note also that the choice [i = aQ is equivalent to an expansion over 
(ln(a 2 <3 2 /A 2 )) -1 rather than over (lnQ 2 /A 2 )~ 1 , i.e. to the change A — > A/a. 

So, let us assume that if one uses a physical scheme with a properly chosen 
subtraction point (i.e., \j? = (k 2 )), then the resultant expansion over (a s /it) has 
coefficients of the order 1 (this is just the situation usually encountered in QED, 
where one has no problems with the momentum-dependence of the coupling constant 
a QED ~ 1/137). If this assumption is valid, then the higher-order corrections may 
be calculated using the following rules: 

1) One may calculate in an arbitrary scheme. The most convenient, in our view, 
is the MS-schemJI^, which is free from spurious terms ln(47r) and je present in 
the minimal-subtraction (MS)-schcmc. The parameter A^jg may be chosen to be a 
fundamental scale of QCD. Note, that A^g is close to Apjj related to a physical 
scheme: Apjj = AtAj^g, where k ~ 2 is almost independent of the vertex chosen to 
define g(k) (cf. Ref.[T5]). 

2) In general, however, there arc no a priori reasons to expect that the 
MS-schemc minimizes the coefficients aik <£> / in Eq. I|2.2[) . If it is known that 

the average off-shellness of lines related to a short-distance subprocess is a 2 Q 2 and 

1 ~ 1 

a»l ( or a < 1), then E® / must be expanded over ln(a 2 Q 2 /K 2 A^jg) rather 



than over 



-i 



[MQ7A^)j 

3) Usually the value of a is not known. However, this value may be estimated 
by requiring that the coefficient aio (or 0,20, if 7o = 0) vanish after the change 
^Ms ~* ^Ajjg/a. 

In this approach all results of the calculations are expressed in terms of the only 
parameter Ao = 2A^jg. However in the expansion (|2.2[) for different processes one 
may use different A^j = Ao/a; (with known a;'s). 

3. Power corrections 

Our derivation of Eq. (jl.ip given in Refs. [SH3 is based on the analysis of Feynman 
diagrams in the a-rcpresentatior|I3 ( se e also Refs. IT8l [19]) . i.e. on the formula 

— H = i da a e m ° ik2 - m ^ (3.1) 

applied to propagators of all lines a of the diagram. After integrations over all fc; 
this gives the representation 



r(Q,p)~ Jl[da a D- 2 {a)G{Q lP ,a) 



x exp 



iQ 2 A{a) + ip 2 I(a) -i V'ay 



to 2 



(3.2) 
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which has many advantages for analysis of the large Q 2 behaviour of T. In particu- 
lar, from Eq. (|3.2p it follows that integration over a region where A(a) > p gives for 
Q 2 — > oo an exponentially damped contribution C[exp(— Q 2 p)\. Hence all contribu- 
tions having a power (0(Q~ N )) behaviour for Q 2 — > oo are due to integration over 
regions where A(a) vanishes. There exist three main possibilities to get A(a) = 0: 

1) short-distance (or small-a) regime, when a ai = a a2 = . . . = a an = for some 
lines 01,02, ...,(T n ; 

2) infrared (or a — > oo) regime, when a ai = a a2 = . . . = a an = oo for a set of 
lines {eri,cr2, ■ ■ ■,cr„}; 

3) pinch regime, when A(a) = for nonzero finite a's. This regime works when 
A(a) may be represented as a difference of two positive terms. 

It is possible also to get A{a) = making up a combination of the three basic 
regimes. In the momentum representation, the first regime corresponds to integra- 
tion over a region k ~ Q, the second one over k ~ p 2 /Q, and the third over k ~ p. 
This means that perturbative QCD is applicable only when the regimes 2,3 and 
the combined regimes either do not contribute at all or give a power suppressed 
contribution compared to that of the pure SD-rcgime. There exists a wide class of 
processes for which the pinch regime does not work (see Ref. [5]), and it is sufficient 
to analyze only the SD- and IR-regimes. In this case it is very useful to visualize 
a diagram as an electric circuit and to treat the parameters a a as the resistances 
of the corresponding lines a. Note, that according to Eq. (|3.2[) . for A(a) = the 
amplitude T lacks its Q 2 -dcpcndcnce. Hence, one must find the subgraphs that 
should possess the following topological properties: when lines of these subgraphs 
are contracted into point (a a = 0) and/or removed from the diagram (a a = oo) 
then the resulting diagram does not depend on Q 2 . Each configuration of this type 
corresponds to some power-behaved C(Q~ JV )-contribution. The power N may be 
easily estimated with the help of the rules ksu ~ Q, ^ir ~ P 2 /Q'- 



responding to the SD-(IR-) integration. Recall that tij — 1 for -0, -0-field.s and the 
curl G^, whereas ti j = for the vector field A^. That is why in QCD (in covariant 
gauges) it is necessary to sum up over external gluon lines of the subgraphs V, S. 

However, for the forward amplitudes (corresponding to inclusive cross sections) 
and for amplitudes of exclusive processes involving colour singlet particles, after 
such a summation the field A^ either disappears (and the gluon lines correspond 
to the curl G M!/ that has a nonzero twist) or enters into covariant derivatives 
Dp = dfj_ — igA^ present in composite operators that naturally arise when the con- 
tribution of the corresponding configuration is written in the coordinate represen- 
tation (cf. Refs.EHZ]). 




(3.3) 




external line of the subgraph V (S) cor- 
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Fig. 2. The configurations for leading contributions to Drell-Yan process with large Q 2 . 

Consider, e.g. the forward amplitude T(r, Q 2 ) corresponding to the total cross 
section of the Drell-Yan process AB — > ji + \x~ X . In this case all the configurations 
responsible for the leading contribution T lead (r, Q 2 ) = 0((Q 2 ) ) have the structure 
shown in Fig. [2^. Here, the subgraph V corresponds to the i?-function, whereas 
subgraphs resulting after contraction of V into point correspond to the function 
/ = /a ® /s ■ The configurations shown in Fig. &jp,c give power suppressed contri- 
butions. Note also that Eq. (|3.3p gives only an upper estimate. This means the 
contribution Fig. [5^ itself apart from the leading contribution, contains also power 
corrections. These corrections appear in the following cases: 

1) If we project the spinor structure of the subgraph V onto 1 rather than 
onto 7^. In the latter case we obtain a composite operator having twist equal to 
2 and in the former one the resulting operator has twist 3. Note that 75,757^ and 
cr^-projections have vanishing matrix elements for spin-averaged amplitudes. 

2) If we expand a bilocal operator M (£, 77) over the local ones, then there appear 
operators 

j times 

|^ (y^i . . . D^) symmctri ^g^..:.g^:^ (3.4) 

having twist 2 + 2j. Each factor g^.^. adds (£ — r/) 2 to the corresponding func- 
tion E(£, 77; 77') and this leads to suppression of the resulting contribution by an 
additional factor 1/Q 2 . 

The expansion of the function E(£, 77, . . .)(£ — r])^ 1 ...(£ — 77)^™ over symmetric- 
traceless structures 77, . . .) [(£ — r]) 2 ] 1 {(£ — r])^ 1 . . . (£ — 77)^'} corresponds in 
the momentum representation to an expansion of the amplitude E(k, Q, . . .) 
over k 2 : 



Eik^.^o + ^ikyEjik,...)^ , 
3=1 



(3.5) 
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i.e. over the off-shcllncss of the particle corresponding to the external line of the 
SD-subgraph V (quark masses are assumed to be zero). In the a-representation 
this corresponds to an expansion of the integrand in Eq. (|3.2[) into a power series 
over Ay = ^ a a . Thus, the coefficient function E for the leading contribution 

corresponds to an on-shcll amplitude. This means that E is formally gauge-invariant 
in each order of PT. However, because of logarithms lnQ 2 /fc 2 present in E(k, Q..), 
taking the limit k 2 — is a rather delicate procedure. Note, that E(Q 2 /fi 2 , . . .) 
in Eq. (jl.ip corresponds to integration over small Ay, i.e. the contribution of the 
region Ay > 1/fJ, 2 must be subtracted off. As a result, one has InQ 2 /^, 2 in place 
of InQ 2 /k 2 , and it is then safe to take k 2 = 0. To maintain gauge invariance, one 
may introduce the cut-off at Ay > l//i 2 using, e.g., the dimensional rcgularization 
d 4 k — > d 4+2E fc(/i 2 )~ e combined with subtraction of poles 1/e. These poles formally 
correspond to ln(/i 2 /fc 2 )| fc 2 =0 (cf. Ref. |2"3")) . 

Thus in the leading power approximation quarks corresponding to external lines 
of the parton subprocess should be treated as the on-shell particles. Real quarks 
are, of course, off-shell. But according to Eqs. (|3.4p . (|3.5[) this phenomenon leads to 
power corrections only. They may be analyzed just in the same way as the leading 
term, although the analysis is more involved. In particular, for each new set of 
operators one must introduce a new function. However some of these functions are 
linearly dependent due to equations of motion. As it was emphasized in a classic 
paper^ 1 , using equation of motion D^^ip = 0, one may get rid of the operators 
containing D^/^ and D^D^ . The resulting operators are built of the fields ip, 'tp and 
covariant derivatives D^. The reduced matrix elements of such operators may be 
identified with the moments of functions which are generalizations of the parton 
distribution functions, e.g. 

<PMK^...JVw...M,>+ fc } \P) = 1 + ( ~ 1} \fa, g (n,k) + (-l) n f- a , g (n,k) 



fag,a(n,k) + (-1)* fa,ga(n, k) + (-1)" fag,a(n, A) + (-!)"+* / s , ffS (n, k) 



1 

fe- 



d.Tidx2dx, 



ii 



fag,a(xi,xr, X 3 )S(xi+X2-X3) + (-l) f a ,ga {x\ J X 2 , X 3 )S(xi-X2-X 3 ) 



+ (-!)*... + (-1) 



n-lk-1 



(3.6) 



where ; denotes covariant differentiation, a denotes the quark flavour and a that of 
antiquark. The function f a g,a (fa,ga) describes a quark (antiquark) with momentum 
X\P and a gluon with momentum XiP in the initial state and quark (antiquark) 
with momentum x 3 P = (x\ +x 2 )P in the final state. The functions f ag . a , fag. a have 
analogous meaning. Such a construction was introduced first in Ref. l25l for operators 

i> {in ■ ■ ■ 9|t fc ) {d^t+i ■ ■ ■ ^Mn+h^n+k+i)} i' useci in tnc analysis of factorization 
in the Feynman gauge. 

Operators of Eq. (|3.6p appear also in configurations 2a if the subgraph V has 
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external lines corresponding to the curl G^ v . Apart from matrix elements (P\0\P) 
these configurations contain also matrix elements (0|G . . . G|0), (0|(?/'V') ■ ■ ■ (ipip)\0) 
(0\G . . . G(ipip)\0), etc. In each order of PT these matrix elements vanish, but 
in QCD, due to nonperturbative effects, these vacuum matrix elements may be 
nonzero. As it was demonstrated in Ref. 1261 these contributions are very important 
for understanding the dynamics of hadrons. The main problem is to generalize the 
methods developed in Ref. 1261 to more complicated amplitudes. 

All the configurations considered above correspond to the SD-regime X(V) ~ 0. 
One must take into account also the configurations 2c corresponding to the com- 
bined SD-IR regime X(V) ~ 0, X(S) — * oo. Physically this regime corresponds to a 
short-distance subprocess accompanied by the exchange of soft quanta between the 
hadrons A and B. According to Eq. (|3.3j) . these contributions also have a power 
behaviour (D(Q~ N ), where N is the number of external lines of subgraph S. If all N 
lines are gluonic, then the corresponding diagrams describe a multipolc interaction 
of hadrons. However, if the quarks are massless, then the subgraph may possess 
quark lines also. The main contribution for the IR-regime is given by the region 
k 2 ~ (p 2 /Q) 2 , where p 2 may be treated as hadronic mass. For massive fields (e.g., 
quarks) the contribution of the IR-regime is damped by the mass present in the 
propagator (fc 2 + m 2 ) -1 if k 2 ~ rn^/Q 2 <C m 2 , i.e. for Q 2 > m A H /m 2 : the IR-regime 
does not work. The mass m in this case works as an infrared cut-off. 

Since the gluons are massless, the IR-regime in PT always works for gluons, 
and there are power corrections due to integration over a — > oo. However, the 
contributions corresponding to the configuration 2c do not factorize in the usual 
sense. This suggests that complete analysis even of the lowest power corrections in 
perturbative QCD is impossible. However, if the exchanged system is coloured, then 
the corresponding contribution should be damped by confinement (nonperturbative) 
effects, i.e. in this case even for m = there exists an IR cut-off M = l/i? C onf ~ 
300-^500 MeV. Thus, if we add to PT a confinement hypothesis, then for a coloured 
system S the contribution of configuration 2c is damped for Q 2 > rn\ 1 /M 2 1 i.e. for 
all hard processes. On the other hand, if the exchanged system is colour-singlet (e.g., 
colour-singlet glue-ball, 7r-meson, p-meson, pomeron etc.), then there are no a priori 
grounds to neglect the configuration 2c. We feel that the methods of the "old" 
hadronic theory, such as the Reggeon calculus and potential models (e.g. the quasi- 
potential approacr|23) will be much more suitable for analysis of these contributions 
than the perturbative QCD methods. Highly instructive in this connection is the 
result of Ref. [28l where it is shown that if one describes the soft exchange by an 
exponentially vanishing quasi-potential, then the soft interactions in initial and final 
states give power (1/pt and l/p?) rather than exponential (exp(— apx)) corrections 
to the amplitude of high-energy wide-angle elastic irp- and pp-scattering. What is 
more, in the available energy range these corrections give an essential contribution. 

Thus, a consistent analysis of power corrections in QCD seems to be a highly 
nontrivial but maybe not a hopeless task. 
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After this paper was essentially completed, we received a preprint by PolitzerJ^ 

where the power corrections are analyzed using the methods similar to ours. There 

is no surprise that the analyzes are similar, because both are based on the clas- 
1041 

sic work!^-* 1 However, there exists also a conceptual difference between the two 
approaches. Our approach is based directly on the analysis of the corresponding 
amplitude T(Q 2 1 p 2 ) in the large-Q 2 region whereas Politzer's approach is based on 
the analysis of small-p 2 behaviour of T(Q 2 ,p 2 ) (i.e. on the analysis of mass singu- 
larities). Both the approaches are (almost) equivalent if one analyzes logarithmic 
m(Q 2 /p 2 )-corrections. However, the power corrections 0(Q~ 2N ) in mass-singularity 
analysis correspond to contributions (p 2 ) N (\np 2 ) k , that simply vanish for p 2 = 0. 
So, we are very sceptic about the main idea of Ref. [29] that a complete analy- 
sis of power corrections may be performed within the mass-singularity approach. 
Incompleteness of this approach reveals itself in the fact that the soft exchanges 
(configuration 2c) are completely ignored in Ref. [29] Another disadvantage of the 
mass-singularity approach is that it incorporates perturbation theory just in the 
region k 2 ~ 0, where one should expect in QCD large nonperturbative effects. In 
particular, within the mass-singularity approach, it seems impossible to understand 
the origin and structure of power corrections to the "e + e~ — > hadrons" process, 
since the related amplitude T(Q 2 ) does not depend (in a massless theory) on small 
momentum variables like p 2 . It should be emphasized that the analysis of these 
effects given in Ref. [55] is the only serious analysis of power corrections in QCD 
and it is based (not by chance) on the operator product expansion, i.e. just on the 
analysis of the large-Q 2 behaviour of the relevant amplitude. 



4. Pion form factor at moderately large Q 2 

As an example illustrating the importance of a detailed study of higher order and 
higher twist effects, let us consider the behaviour of pion electromagnetic form factor 
F n (Q) for moderately large Q 2 

During the last 3 years a definite progress has been made in understa nding of 
the asymptotical Q 2 — ► 00 behaviour of F^^Q) in the QCD framework EHHEDHHH] 
The main result here is the proof that in a region where the power corrections may 
be neglected the form factor may be written in a factorized forrrP^ - ^] 
1 1 

F *(Q) =^2 J dx J d 3/V*(j/)M 2 iMfl)a«(MH)) 


x E(Q 2 /fi 2 , Vx/fJ 2 , x, y, a a (fi R )) ip(x, /i 2 , fj.% a s (^ R )) , (4.1) 

where tp{x) is the wave function describing the splitting of the pion into a qq-st&te 
and E/Q 2 is the amplitude of the short-distance subprocess j*qq — > q'q'. Note, 
that, in general, the normalization parameter fi R of the R-operation for ordinary 
UV -divergences may differ from the splitting parameter /1 that separates small and 
large momenta. The latter may be treated also as the normalization parameter for 
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composite operators. The moments of the function (p(x) are equal to the reduced 
matrix elements of the twist-2 operators ip r y^D n ip. The SD-amplitude, as usual, is 
a series expansion over a s (fin) (see Fig. [3]): 

vnirfli 2 2/2 , \\ 2Tra s (n R )C F xQ 2 

E(Q /n , V R l\i , x, y, a„(ji R )) = — ( x Q 2 )( xy Q2) i 1 + °( a °)} > 

(4.2) 

where Cp = 4/3, N c = 3. The factor xQ 2 in the numerator of of Eq. (|4.2[) is 
due to the trace over Dirac 7-matrices. Note that E(x, y) is rather singular for 
x,y —> 0. Hence, the main contribution is given by integration over small x,y. 
In the next order the most singular terms are hi{xyQ 2 / y 2 R ) and ^(lna;)(ln y), 
(In Q 2 1 fi 2 ) ■ (In xy) / (xy) . The first term is given by divergent parts (notice ^lr) and 
the others by convergent ones. Thus, to minimize the a s -corrections we must take 
fj? ~ xQ 2 = yQ 2 and [i 2 R ~ xyQ 2 , where x (y) is the average value of x (or y) : 

If f(x) ~ 5 (x — 5) (noninteracting quarks), then x — ^. However, for very broad 
functions, e.g. for tp(x) ~ [(x(l — x)] R with fl < 1, we have a very small value 
x ~ exp(— l/R). 

In the problem investigated we en- 
counter just the same mass scales 
(m T ,m g and M = l/R con i) as in deep 
inelastic scattering. So, it seems natu- 
ral to expect that Eqs. (|4~Tj) . f4~2|) . for 
the appropriately chosen wave function 
tp(x) must provide a good approxima- 
tion for F 7r (Q) in the region Q > 1 GeV. 

The wave functions ip(x, y 2 ) are in gen- Fi s- 3 - Dia S ram for asymptotic term of pion 

. form factor. 

eral unknown. In perturbative QCD one 

may calculate only their evolution with growing fi 2 . In particular ip(x,y 2 ) — > 
6/ w x(l — a;) as /i 2 — * 00, where / ff = 133 MeV. Presence of the /^-factor is due 
to the normalization condition 

1 

cp(x,(j?)dx = U. (4.4) 


It is clear, however, that for /i < 1 GeV the wave function <p(x, /1 2 ) may strongly 
differ from its limiting form. For noninteracting particles <p(x) ~ 5(x — 1/2). When 
the interactions are switched on, the wave function broadens. The width T of ip(x) 
may be estimated as T ~ (Ei nt /m q ) 2 . Hence, for heavy mesons (e.g. for J/ip or Y- 
particles) <p(x) is rather narrow, since Ei nt ~ M = 300 -f- 500 MeV and m g > 1 GeV. 
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On the other hand, for pions the wave function must be very broad, because 
m u ~ 4MeV, ma ~ 7MeV, i.e. pion must be treated as an ultra-relativistic system. 
To obtain a more accurate estimate of the width r for such a system we assume 
that the (soft) Bethe-Salpeter wave function xpfti, fe) is exponentially damped for 
moderately large spacelike k 2 (i = 1,2): 

X p(kuk 2 )~-±- s e X p(k!/M 2 ); —kf > M 2 . (4.5) 

K 1 K 2 

The exponential damping is suggested by the observed spectra of particles produced 
in high-energy hadronic reactions. F or our wave function tp(x) (which may be ob- 
tained from Xp(^1j^2) by integratiorJ^ESl over fc_ = fc — ^3 an d k±), the choice 
(|4~5|) gives 

{exp[-ml/xM 2 } ; x < 1 

(4.6) 
exp[-ml/(l - x)M 2 ] ; (1 - x) < 1 . 

Thus, <p(x) is very close to f n everywhere outside the regions < x < rn^/M 2 
~ 10~ 4 ; < 1 — x < m 2 d /AI 2 ~ 10~ 3 . In these regions f{x) vanishes rapidly. Note 
that for such a wave function x ~ 10 -3 -j- 10 -4 , i.e. the main contribution into F n is 
given by the region where the gluon has a catastrophically small off-shellness x 2 Q 2 , 
which for Q 2 < 100 GeV 2 is much smaller than the value \k 2 \ ~ 0.1^-0.3 GeV 2 , where 
the confinement effects must be taken into account. Thus, for a broad wave function 
short distances do not contribute, in fact, and Eq. ()4.2|) is unreliable. In particular, 
it is not justifiable to neglect the power corrections that may really have a form 
(M 2 /(fc 2 )) ~ (M 2 /xQ 2 ) rather than simply M 2 /Q 2 . We assume that the confine- 
ment effects eventually remove the infrared singularity from the "hard" quark and 
gluon propagators 1/xQ 2 and 1/xyQ 2 . So, we change 1/xyQ 2 — ► l/(xyQ 2 + 2M 2 ) 
- 1/(^Q 2 + ((Ar x - fc^) 2 )) and 1/xQ 2 -» l/(xQ 2 + A/ 2 ) - 1/(.tQ 2 + (fc 2 ^)). The 
connection between M 2 and (kj_) is a pure mnemonics and must not be understood 
too literally. However, as an order-of-magnitude estimate this connection must be 
true. So, we should expect that M 2 ~ 0.1 -r 0.3 GeV 2 . 
As a result, we have in place of Eqs. (|4.1|) . (|4.2j) 

l 

?AArn\ _ 2rrCF (a a i \ i \ a s {n 2 R )xQ 2 



^ A (Q) = -T7 il / dxdy<p(x)cp(y) 



N c J ' rya/ (xQ 2 + M 2 )(xyQ 2 + 2M 2 ) 

o 

x{l + 0(a s )}, (4.7) 

where AA stands for projection onto the "axial" operators tpj^j^D""^. 

From Eq. (|4.7| it is clear that for not too large Q 2 the contribution of the soft 
region x ~ is damped at x ~ Al 2 /Q 2 , whereas the wave function damps only the 
region x < m 2 /M 2 . Thus, up to Q 2 ~ {M 2 /m q ) 2 > 10 3 GeV 2 the magnitude of the 
pion form factor is determined by the value of M 2 , i. e. by the confinement radius. 
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The main contribution into the integral in Eq. (|4.7|) is given by the region 
xyQ 2 ~ 2M 2 . To minimize the a s -correction we should take [i 2 R equal to the average 
off-shellness of the gluon: fx R = xyQ 2 + 2M 2 ~ 4Af 2 , 

47T 

« 8 (4)^ 8 (4M^ 9h(4M2/A2pH) . (4.8) 

It should be realized that Eqs . (|4.7p , (|4.8[) are meaningful only if a s (4M 2 ) /V < 1 , 
i.e. for AM 2 /A 2 > 50 (in this case a s (4M 2 )/7r ~ 0.1). If A/ 2 - 0.2 GeV 2 , then Eqs. 
(|477|) . P~g]) may work only for A PH ~ 100 MeV. Note that this is just the value 
preferred in Ref. Hfo] We emphasize that the authors of Ref. just have taken into 
account power corrections. In standard analyzes of deep inelastic data (neglecting 
higher-twist effects) larger values of A are usually obtained. It is known, however, 
that if one includes in a phenomenological analysis the effects of higher twists, then 
it is possibkJ^to describe the data using an arbitrarily small A. 

Apart from power corrections related to the primordial transverse momentum 
of quarks (which correspond to operators involving the curl G^ u ), there exist also 
power corrections due to twist-3 operators ^^D n ip and ip^^a^ v D n ip. In the large- 
Q 2 limit their contribution has an additional factor A 2 /Q 2 compared to the contri- 
bution of the twist-2 operators V'TsO '^D n 't/j . Note, however, that A is anomalously 
large 

- TR 2 

(0\d l5 u\P) = %U — - if, ■ (1-8 GeV) , (4.9) 

m u + rrid 

i.e. for Q 2 < 6 GeV 2 these operators cannot be neglected^. For the pseudoscalar 
, 075D rl '0-operator we have 



P p 47ra s (4M 2 ) . 
tt (Q) = TT °f / <PP{x)<pp{y)dxdy 



A. 



1-x 



{xQ 2 + M 2 )(xyQ 2 + 2M 2 ) 



{l + 0{a s )}, (4.10) 



where ipp{x) = \ip(x) ~ A/^. 

It should be noted that for M 2 = the amplitude E^ pp \x, y) is as singular at 
x ~ as 1/.t 2 . As a result, integration over x gives an additional factor Q 2 /M 2 that 
compensates the absence of the Q 2 -factor in the numerator of Eq. (|4. 1 0|) . In other 
words, in the region Q 2 <C (M 2 /m q ) 2 the contribution F PP (Q) has l/Q 2 -behaviour 
rather than 1/Q*. Moreover, F PP (Q) has an additional large factor (A/M) 2 > 10 
compared to F^ A . The same factor has also the -F!j p -contribution ((TP) stands for 
tyl&tivty ® ?/'75V'-P ro j ect i on ) : 



h (Comm. 2009) B. V. Geshkenbcin and M. V. Terentev, in Phys. Lett. B 117, 243 (1982), showed 
that one should also include tp-fs^y^D j_D n rp operators. 
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Fig. 4. Dependence of pion form factor predictions on the confinement parameter M. 



Fit (Q) 



2na s {4:M 2 )C F 



1 - 



ip T (x) (fp(y)dx dy 
(xQ 2 + M 2 )(xyQ 2 + 2M 2 ) 

{l + 0(a s )} 



(l + x)Q 2 _ y(l + x)Q 2 
xyQ 2 + 2M 2 



(4.11) 




xQ 2 + M 2 

Note that for M 2 = the am- 
plitude E TP has the l/x 3 -singularity. 
However, using the equations of motion 
it is possible to show that the func- 
tion ipt(x) has an additional x- factor 
for x ~ 0. In particular, if tpp(x) ~ 
A/tt, then (Pt(x) ~ Xf v x(l — x). As a 
result, Fj p (Q) ~ 1/Q 2 in the region 
Q 2 ^ (M 2 /m q ) 2 . This contribution is 
negative and small for Q 2 > 4GeV 2 . 
The curves for F n given by the sum of 

Eqs. P~7|) . (gUnj), (|SIIII> have a right form (Fig. gj) and for M 2 
A = 100 MeV they are close to existing experimental data. 

It is easy to realize that since the passive quark in our case has a very small 
fraction of the pion mom entum ( "wee" parton) , we deal really with the mechanism 
proposed by FeynmarPSI to explain the power-law fall-off of hadronic form factors. 
Thus, we might have considered th e diag ram shown in Fig. [5] and write F^{Q) in 
the standard bound state formalisrrJ23311 



Fig. 5. Pion form factor in light-front formalism. 



0.1 ^-0.2GeV 2 , 



F„(Q) 



d.r 



x(l — x) 



d 2 k± 0(1 — x, k±) 0(1 — x, k± + xq±) 



(4.12) 



Note that according to the Bcthe-Salpetcr equation the diagrams shown in Figs. 
3 and 5 are equivalent up to O(a s )-corrections. 
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If the function 4>(x, fcj_) is that given by Eq. (|4.6p then, performing k± -integration 
in Eq. (|4.12|) we obtain 



FAQ) ~ J da; exp 



xQ 2 



2Af 2 (l-:c) (1 - x)Af 2 



(4.13) 



whence it follows that for Q < M /m 2 the main contribution is given by the region 
x ~ M 2 /Q 2 . If the function <^(a;) (i.e. fex) integrated over k±) behaves like x R 
for x ~ 0, then F 7T (Q) ~ (Q 2 ) -1- ^. Our choice l|4.5[) corresponds to i? ~ 0, and as 
a result F 7r (Q) ~ l/Q 2 . If one assumes that 

X p(fci,...,fc„) ~cxp{^fc 2 /A/ 2 } (4.14) 

for a system composed by n valence quarks, then y(x) ~ x n ~ 2 for a; ~ and thus 

F (n) (Q)~(Q a ) 1_B . (4-15) 

This r elation corresponds formally to the well-known quark counting rule 
our case, however, this rule has nothing in common with short 
distances and scale invariance. The short-distance mechanism proposed in Ref. [40] 
to explain QCR, according to our estimate, works only for Q 2 > 10 3 GeV 2 . In an 
intermediate region one must take into account the fact that the contribution of the 
Feynman mechanism is damped by the Sudakov form factor of the active quark. 
Thus , on e must multiply the curves shown in Fig. [4] by the Sudakov QCD form 



factoiSH 



S(Q 2 ,M 2 ) = exp 



2CV 
' b 



2 3\ / lnQ 2 /A 2 \ , Q 2 ^ 



(4.16) 



For A = 100 MeV, M 2 = 0.22 GeV 2 this gives the curve shown in F ig. RH In the 
region Q 2 = l-j-4 GeV 2 there is good agreement with experimental data!^ Decrease 
of Q 2 Fk(Q 2 ) for Q 2 > 10 GeV 2 is due to the Sudakov form factor. In the region 
Q 2 > 100 GeV 2 the short-distance regime begins to work. In this region the average 
off-shellness of the gluon increases, fj, 2 grows, and as a result the wave function 
becomes narrower: 



<p(x,fi^) ~ ( x (l-x))^r ln ^~ lnln ^\ , (4.17) 

and this, in turn, damps the contribution of the Feynman regime. For Q 2 > 
10 3 GeV 2 one may neglect the F^ PP ^ contribution and use Eqs. ([4.1ft . (|4.2[) with the 
wave functions <p(x) ~ (x(l — x)) . The asymptotic formula tp{x) — &f w x(l—x) 

may be used only for Q 2 > 10 20 GeV 2 . 

Thus, in the region Q 2 > 10 3 GeV 2 begins the asymptopia, and F 7T (Q) is again 
given by the quark counting rules. In this case they are due to the short-distance 
scale invariance, as expected in the pioneering works by Matveev, Muradyan, and 
TavkhelidzJS^ and Brodsky and FarrarP^ 
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Q 2 F n (Q*) (GeV 2 ) 




1 4 10 100 1000 Q 2 (GeV 2 ) 



Fig. 6. Prediction for pion form factor with Sudakov form factor included. 

Quark counting rules for ultra-relativistic systems were considered first by 
Terentyev.SS However, he used constituent quark masses, m q ~ M . In this case 
the range of applicability of our analysis (M 2 <C Q 2 <C M 4 /m 2 ) is zero. We insist 
on using the current quark masses (m q ~ 4 — 7McV) in the wave function (|4.6|) . 
Note that the eventual IR cut-off in our analysis is of an order of M (i.e. of an order 
of constituent quark mass) in accordance with the common wisdom. 

The last but not least observation is that the magnetic proton form factor in 
QCD is negative for narrow wave functions^! like ip(x\ ,X2,x%) ~ ^[5(^ — 1 /3) and 
positive for the broad ones, e.g. for that given by Eq. (|4.f 4[) . 

Summarizing this section, we may conclude that although our analysis is semi- 
phenomenological and some assumptions are very crude, it is clear, nevertheless, 
that a consistent treatment of power "corrections" (in fact, they give the main 
effect) is the main problem for perturbative QCD of hard elastic processes in the 
now (and, perhaps, forever) available energy range. 

5. Soft processes and perturbative QCD 

The main fraction of the total cross-section at high energies is due to the processes 
with small transverse momentum (soft processes). These are the elastic and quasi- 
clastic processes in the diffraction region ( C 1 (GcV/c) 2 and multiple production 
processes with low pt- pt <C 1 GeV/c. The conventional phenomenology of processes 
in this region is the Regge-Mueller picture. Till now, however, it was not clear 
whether perturbative QCD can give any information about these processes. Below 
we discuss this problem analyzing a process 12 — > 1'2' in the region s ^> t, ?7i 2 lac i r - 
We will assume also that the i-channcl is flavour nonsinglet. This assumption greatly 
simplifies the analysis. 

For scalar gluons, i.e. in a Yukawa-type field theory soft processes have been 
studied fO years ago SSI it was shown that summation of all logarithmic terms 
(logS) N coming from the short-distance integration (regime f, see Section[2]) gives 
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the following representation 




(5.1) 



F ± (S,t) = - 




r(j + 1) sin(Trj) 



\S\ j (e inj ± 1) 



(5.2) 



where ± stands for signature, C, v and B are some matrices (to be discussed below); 
e.g., B = B a b, and a,b~S, V, T, A, P arc structures appearing in the Fierz identity 
applied to factorize the spinor structure of the relevant contributions. According 
to this representation (Eq. (|5.1[) ) the Mellin transform /^(j, i) possesses moving 
(t-dependent) Regge poles due to zeros of Det[l — B(j,t)v(j)]. It has also fixed 
(t-independent) singularities in the complex j-plane accumulated in the function 
v(j). The type of fixed singularities depends on the ultraviolet asymptotics of the 
effective coupling constant. In particular, in a fixed point theory (where g(fJ.) — * go 
as /i — > oo) the function v(j) has square-root branch points, the position of which 
depends on go, i.e. on the asymptotical value of g(fi). On the other hand, in an 
asymptotically free field theory v(j) has the infinite number of poles condensing to 



Let us now discuss briefly the derivation of Eq. (5.1). Consider a particular 
diagram of a binary process 12 — > 1'2' in the region s ^> \t\, mhadr The Mellin 
transform of its contribution has the following structure in the a-representation 



/ Y[da a D- 2 (a)g(j,t,a)\A(aW[e(A)±e(-A)}cxp[iJ(a,t,m 2 )}, (5.3) 



where g(j,t,a) is a polynomial in j (it corresponds to the function G in Eq. (|3 . 2[) ) 
and A is the coefficient in front of the large variable S = s — u. As is well known, 
the asymptotical behaviour of F(S,i) for large S is determined by the rightmost 
singularities of its Mellin transform f(j,i). These are poles j — N generated by 
integrations corresponding to the regimes l)-3) discussed in Section [3] However, 
using Eq. (|3.3jl it can be shown that in Yukawa theory the IR-regime gives only 
non-leading poles at j = —1,-2,.... Furthermore, the pinch regime contributes 
only to the negative signature amplitude F~ (S, t). Thus, for F + (S, t) it is sufficient 
to consider only the poles due to the short-distance regime. 

According to Eq. (|3.3p . the leading poles (at j = 0) are due to the subgraphs 
Vi with 4 external lines. We recall that V* should possess the property that if it 
is contracted into point, the diagram becomes S'-indcpendent (i.e. Vi must be an 
S'-subgraph). The most general configuration is shown in Fig. [71 Note that in general 
the SD-subgraphs Vi may be 2-particle-reducible, i.e. they may contain smaller S- 
subgraphs with 4 external lines and the total singularity due to the SD-regime of 
Vi may be a multiple pole j~ Ni - It makes sense to treat a particular diagram as 



j = o. 



(Eq. (32)): 



<7 
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Fig. 7. Illustration of Equation 15. II) . 

a ladder composed by 2-particle irreducible blocks kj. Then the maximal value of 
Ni is determined by the number of kj 's inside Vi (and also by the number of the 
UV-divergent subgraphs inside Vi). The contribution fv(j) of each S'-subgraph V 
may be represented as a sum of two terms fy = f v ° le + fy S - The first term (fy 0l °) 
is due to integration over the region a a = X v < 1 //i 2 and the second one is due 
to that over the region Ay > l//i 2 . This procedure corresponds to a subtraction of 
the pole due to the small Ay integration. However, if V is composed by two or more 
fcj's then / rcg may also possess the poles at j = due to the SD-integration for a 
smaller subgraph Vi C V. Thus, one must represent fy S as fy s V^ po ° + g V* eg , and 
so on. An example of such a decomposition is shown in Fig. [8l where the pole parts 
are circled by the thin (red) line and the regular ones marked by the slashed (blue) 
lines. Note that Fig. [7] is really a decomposition of the whole diagram. Summing 
over all diagrams we obtain (in the coordinate representation): 



{2nf8 i (p 1 + P2 -p' 1 -p' 2 )F(S 1 t)^Y, j CWi;zi,yi) (5.4) 

n-l n=l 

, yi ; Zi , Wi) B (zi , Wi ; x i+1 , y i+1 ) } v (x n , y n ; z n , w n ) C(px , p[ ; z n , w n ) [d%„] 

i=l 

where [dx]i = dxidyidzidwi (see Fig. [7]), and the functions C, v, B are given by the 
following matrix elements 




Fig. 8. The leading terms of a particular box diagram. 
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C a (pi,j/ 1 ;xi,yi) = (p[\S + T(: ip(x 1 )T a %p(y 1 ) : S)\p[) , (5.5) 
B a b(zi, w l ;x l+ i,y l+ i)^R fl 2{0\S + T(:rp(z i )r a ^(wi): :$(xi + i)Tbip(yi+i): S)\0) ,(5.6) 
v(x t ,yi; w ?+ i)=P Al 2 (0\S + T(:fj(x i )T a r](y i ): :f)(z i+1 )T b r)(wi +1 ): S)\0) , (5.7) 
T a , r fc = 1, 7 P , 07^, 75, 757^ . (5.8) 

Here, : : denotes the usual normal product; 77, fj are the spinor currents (e.g. 
7] = S + (SS/S'4'))) S is the S-matrix; P^i means that Ay < 1/p 2 for each diagram 
V contributing to v, and R^i means that Ay > l/^ 2 f° r au leading ^-subgraphs 
having lines related to B or C. 

If we expand B or C into the Taylor series over £ = x — y, £ = z — w, 

00 

C a ( Pl ,p[;x 1 ,y 1 ) = exp(irX 1 )R^-(p[\S+T(O a ^... H S)\p[)^- ■ (5.9) 

3=0 J ' 

B a b(z ii w i ;xi +1 ,yi +1 ) = (5.10) 

OO OO 1 

^EE^m ' -CClti- ••Cr+ fc i(0|5+T[O aI/1 ...^(Z i )O i , Ml ...^(X i+1 )5]|0), 

(where r = pi — p^, X t = z< + y< ; £ 4 = Zi +" l ) 1 then i? M 2 provides the renormal- 
ization recipe for the resulting composite operators O aVl ... v . = ^T a d Vl ■ ■ -d v .ip. It 
should be emphasized, however, that in addition, subtracts from B, C also the 
contributions due to integration over small Ay-parameters for S-subgraphs that do 
not contain the vertices ipTd n ip corresponding to composite operators. 

Just like for hard processes, only the lowest-twist operators (i e. the traceless- 
symmetric part 0{ avi ... v .y of O) give the leading contributions. Using the transla- 
tion invariance of the functions (|5.5(> — ()5.7|) . integrating over £j, Xi and Zi and 
summing over n (n is the number of SD-intcgrations), we obtain 

F(S,t) = J2 l -^-C(j,t) ■ [1 - v^B^t)]- 1 v(j) ■ C(j,t) (5.11) 

where t = r 2 , S = 2{PQ), P = Pl + pi, Q = p 2 + p'i- The functions B(j, t), C(j, t) 
and v(J) are given by 

R^(p[\S + T(O a ^... H SM)=C(jMP a P» 1 ---P f > j } + 0(r fi ) (5.12) 
J dXe irX R^(0\S + T[O a „ 1 ...„ J (Z i )O b ^--' 1 '°(X i+1 )S}\0) = 

J d(X - Z)d^dC exp[zr(X - Z)] r2 = v(X - Z, f , ■ ■ ■ ' ' ■ 0* = 

iKiK^ 1 1 ---^ } } + 0(r,) J (5.14) 

where w(X — Z, £, £) is defined by 

v(x,y,z,w)=v(x + £,X-£,Z+£,Z-± \ =v{X-Z,t,0 (5.15) 
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and 0{r v ) denotes terms containing r v . These give zero contribution into Eq. (|5.1ip 
because (Pr) = (Qr) = 0. Note that Eqs. (|5.1[) . (|5.2[) give just a Mellin transformed 
version of Eq. (|5.1ip . 

To construct the functions B, C and v one must apply first the R-opcration for 
ordinary divergent subgraphs (this procedure is characterized by the renormaliza- 
tion parameter /ir) and then the operations P^ and R^2 = 1 — P^i that separate 
small and large A-parameters (this corresponds to splitting of the mass logarithms 
hiS/p 2 into "short-distance" (InS/fj, 2 ) and "long-distance" (ln/x 2 /p 2 ) parts). The 
whole amplitude F(S,t), of course, must be independent both of fi and \xr. The 
/Xfl-independence of F leads to a standard renormalization group equation 

(fiRd/duR + (3(g)d/dg - 4^)v(j, fi R , g, fx) = . (5.16) 

The subtraction procedure for the problem considered is more complicated 
than that for hard processes. This is mainly due to the fact that, for soft processes, 
we deal in general with the configuration (Fig. [7]) that has several non-overlapping 
SD-subgraphs V\ , . . . V n . We recall that for hard processes we always have a con- 
figuration with the only SD-subgraph (see, e.g., Fig[2k). Straightforward analysis 
gives the following equation for C 

= C' = j(j)C + (1 - Bv)~\l - Bv)' , (5.17) 

d/x 

where all functions entering into Eq. (|5.17p : B, C and v depend on \i. The sec- 
ond term in Eq. (|5 . 1 7[) is just due to the additional subtraction discussed above. 
The function is the ordinary anomalous dimension of the composite opera- 
tor O avi ... v .. In our case it is convenient to single out from 7 the terms singu- 
lar at j = 0. It can be shown that these terms are proportional to v, i.e. that 
= c(j) + b(j)v(j), where c(j) is regular at j = 0. 
The equation similar to Eq. (|5.17|) can be obtained also for (1 — Bv): 

- (1 - Bv)' = [2 7 (j) + bv(l - Bv) + Bv'][l - Bv}. (5.18) 

Requiring that $ pole (j), the sum of the leading poles, does not depend on \x we 
obtain that <I> polc (j), must be regular at j = 0: 

M A$Pole (i) = r(i) (5 . 19) 

d/Lt 

where r(j) is some function regular at j = 0. Using Eqs. (|5. 17115. 19|) we obtain the 
equation for v: 

v' + 2jv + bv 2 = -r. (5.20) 

The meaning of r(j) becomes clear, it is just the residue of v(J) at j = 0, because 
v' = jv. 

It should be remarked that Eq. (|5.20[) differs from its analogue given in Ref. H51 
because of another choice of the i? 2 r operation. 
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Using Eq. (|5.20|) one can sum up all the poles at j = due to the SD-regime 
of all possible S-subgraphs (i.e. to sum all log N (S/p 2 ) contribution). The solution 
of Eq. (|5.19p has square root branch points in the complex j-planJS51 (see also 
Ref. [46]) . However, v(j) has also poles due to divergent subgraphs. These poles 
(i.e. log JV (S'/^|,)-contributions) are summed by Eq. (I5.16|) . If we take [i = and 
combine Eqs. (|5. 16|) . (|5.20[) . we obtain 

dv 

P(g)-^ = (j - 2 7 - 4 7 v> " bv 2 - r , (5.21) 

where v, 7 and r depend on j, and {—j) is the canonical dimension of v. In the 
lowest order of PT b = 1, r ~ 7 ~ 7^ ~ g 2 , (3(g) ~ g 3 , and the solution of Eq. 
(|5.2ip has condensing poles at j = ) 45 | 46 | 

Summarizing the preceding discussion, we conclude that if one assumes that the 
asymptotical behaviour of F(S, t) is given by the sum of the leading terms of all 
contributing Feynman diagrams, then F(S, t) has for large S a Reggc-type behaviour 
F(S,t) ~ C 2 (t)S a ^ since its Mcllin transform has just a i-dependent singularity 
at j = a(t). To find the function a(t) explicitly, one must solve the equation 

Det[l - B(j, t, fj, R , g, fi, m)v(j, (j,r, g, fi)] = . (5.22) 

It can be shown that Eq. (|5.18[) guarantees that a(t) does not depend on fi and fin: 

a(t) = ^m 2 q /t, t/^ 2 ,g(^ 2 )) = ^{m 2 Jt, 1, g(t)) (5.23) 

Hence, one may try to calculate the Regge trajectories in the region where g(t) is 
small, e.g. in QED, where a QED « 1/137, or in QCD for sufficiently large t. However, 
there arises a question whether Eq. (|5.1ip is valid in vector theories. 

In QCD one encounters the complication discussed in Section |3l First, a SD- 
subgraph Vi may have an arbitrary number of external gluon lines. But if the 
t-channel is colour singlet, then the only change is 

$(x)Tfl>(y) -» ^(x)rPexp [ ig [ A^z)dz^y) I i>{y) (5.24) 



for all bilocal operators entering into B- and C-functions. For local operators this 
corresponds to the change — > = — igA^. The second complication is due 
to the IR- regime (soft exchanges, see Fig. [9]). However, just like for hard processes, 
if the ^-channel is colour singlet, then the sum of all soft exchanges give only power 
corrections in each order of perturbation theory. Thus, all terms responsible for 
the leading power contribution have the structure of Fig. [7] and as a result, we 
get Eq. (|5.1ip . In other words, if we sum the leading j singularities of all relevant 
Feynman diagrams in QCD, we obtain a Regge-type picture for binary processes 
and, hence, a multiperipheral picture for the multiple production at low pt- 
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We recall, however, that we have 
discussed above only the flavour nonsin- 
glet, positive signature amplitude F^ s . 
For F^s the pinch regime (see Sec- 
tion [3]) also gives leading j-poles for 
non-planar diagrams. It is known, how- 
ever, that the non-planar diagrams have 
an additional colour factor (l/iV c ) 2 = 
(1/3) 2 . This suggests that the pinch 
contributions in QCD must be sup- Fi §- 9 - The samc as Fi S- El but for the case of 
pressed. There exists also an experi- ^ D ' 

mental evidence in favour of this suppression: the well-known signature degeneracy 
of the Regge trajectories. 

For flavour singlet amplitudes F$ ("vacuum" exchange) the poles generated by 
the pinch regime are at j = 1 rather than at j = due to the 2-gluon intermediate 
states and after summation one obtains for Fa a square-root branch pointEl at 
7 = 1 + 0(g 2 ). This suggests that the pinch regime plays a highly important role in 
formation of the Pomeron singularity. 

The most intriguing possibility is to utilize the asymptotic freedom of QCD for 
a calculation of the Regge trajectories and of the resonance masses in the region of 
large t (see Eq. (|5.23p ). Note, however, that the function B(t) describes the long- 
distance dynamics, i.e. by its construction, B has an UV cut-off but there is no IR 
cut-off. This means that if the IR region of integration gives a sizeable contribution, 
one must (in some way) take into account nonperturbative effects. It seems that 
the most effective tool here is the method proposed in Ref. HSj This and related 
problems are under investigation now|f| 
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